Multi-particle States in Spin-1 Chain System CsNiCls 



M. Kenzelmann\ R. A. Cowley\ W. J. L. Buyers^'^, R. Coldea''■^ J. S. Gardner^, M. Enderle^ D. F. McMorrow'', 

S. M. Bennington^ 

(1) Oxford Physics, Clarendon Laboratory, Oxford 0X1 3PU, UK (2) Neutron Program for Materials Research, National 
Research Council of Canada, Chalk River, Ontario, Canada KOJ IJO (3) Canadian Institute for Advanced Research (4) Oak 

Ridge National Laboratory, Solid State Division, Oak Ridge, TN 37831, USA (5) ISIS Facility, Rutherford Appleton 
Laboratory, Oxon 0X11 OQX, UK (6) Technische Physik, Ceb. 38, Universitdt des Saarlandes, 66123 Saarbriicken, Cermany 
(7) Condensed Matter Physics and Chemistry Department, Ris0 National Laboratory, DK-4000, Roskilde, Denmark 

(February 1, 2008) 

A continuum of magnetic states has been observed by neutron scattering from the spin-1 chain 
compound CsNiCla in its disordered gapped one-dimensional phase. Results using both triple- 
axis and time-of-flight spectrometers show that around the antiferromagnetic point Qc ~ tt, the 
continuum lies higher in energy than the Haldane gapped excitations. At 6 K the integrated intensity 
of the continuum is about 12(2)% of the total spectral weight. This result is considerably larger 
than the 1-3% weight predicted by the non-linear sigma model for the 3-particle continuum. 



PACS numbers; 75.25.+Z, 75.10.Jm, 75.40. Gb 

The excitations of one-dimensional (ID) Heisenberg 
antiferromagnets have attracted much experimental and 
theoretical attention ever since Haldane |l[ predicted that 
the excitations of integer spin and half-integer spin chains 
are different. Half- integer spin chains have no spin gap 
and exhibit a spinon continuum extending to zero en- 
ergy I^J^ . For integer-spin chains the spectrum is at low 
temperatures dominated by well-defined single particle 
excitations corresponding to a triplet of spin-1 particles 
that exhibit a large energy gap |^,^. These have been 
thoroughly studied by neutron scattering from spin-1 
chain compounds and it has been shown that the well- 
defined excitations largely exhaust the total scattering 
1^. Nonetheless, a continuum of multi-particle scatter- 
ing was predicted by Haldane and others Q, but as yet 
no experiment has established its existence. To search for 
such a continuum we have measured the inelastic neutron 
scattering of CsNiCla in its paramagnetic ID phase us- 
ing time-of-flight and triple-axis spectrometers. We find a 
neutron scattering continuum at energies higher than the 
well-defined single particle excitation, whose integrated 
intensity is bigger than predicted for the multi-particle 
scattering. 

CsNiCls is magnetically a quasi ID spin-1 chain com- 
pound with the Ni chains along the c-axis of the hexag- 
onal unit cell. Each unit cell contains two Ni atoms and 
the Hamiltonian is 



chain 



plane 



= J ^ 5i • + J' ^ 5i • 5j ~Dj2iSn' ■ (1) 

i <i,j> i 

The exchange interaction J — 2.28 meV along the c- 
axis is much stronger than the exchange interaction in 
the basal plane J' = 0.044 meV The weak Ising 

anisotropy D = A fieV is small enough that CsNiCls is 
a good example of an isotropic Heisenberg antiferromag- 
net. Below Tn = 4.85 K the interchain coupling causes 
long-range ordering of the magnetic moments. Above Tn 



CsNiCla is in a ID magnetic phase for which the mag- 
netic exchange interaction along the c-axis dominates the 
spin dynamics, except close to Tn- 

The sample of CsNiCls was a single crystal 20 mm x 
5 mm X 5 mm and was mounted with its (hhl) crystal- 
lographic plane in the principal scattering plane. Inelas- 
tic neutron scattering experiments were carried out with 
the time-of-flight spectrometer, MARI, at the ISIS fa- 
cility of the Rutherford Appleton Laboratory and with 
reactor-based triple-axis spectrometers, DUALSPEC at 
the Chalk River Laboratories and RITA at Ris0 National 
Laboratory. 

For the experiment using MARI the incoming neutron 
energy was 20 and 30 meV. The energy resolution was 
0.35 and 0.4 meV, as determined from the full width at 
half maximum (FWHM) of the quasi-elastic peak. The 
resolution in wave- vector transfer at zero energy transfer 
was typically 0.02 A^^ along the c'^-axis and along the 
[110] direction and up to 0.19 perpendicular to the 
scattering plane if only the central detector bank was 
used. Both the energy and the wave-vector resolution 
improved with increasing energy transfer. The measure- 
ments were performed with the sample at temperatures 
between 6.2 and 12 K. 

The measurements at the DUALSPEC triple-axis spec- 
trometer were performed at 8.5 K with a fixed scattered 
neutron energy of 14.51 meV and with a graphite fil- 
ter to absorb the higher order reflections from the py- 
rolytic graphite monochromator and analyzer. The coUi- 
mation from reactor to detector was 0.65°-0.6°-1.4°-2.0° 
and gave an energy resolution of 1.94 meV, as deter- 
mined from the quasi-elastic peak. In the case of the 
RITA triple-axis spectrometer, supermirror guides and 
a rotating velocity selector produced a variable energy 
monochromatic neutron beam with suppressed higher or- 
der contamination. The scattered neutrons were ana- 
lyzed using the SoUer geometry with a fixed energy of 
5 meV and a cooled beryllium filter between the sam- 
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pie and analyzer. The effective collimation was 1.45°- 
1.45°-l°-2.0° and the energy resolution was 0.35 meV as 
determined from the quasi-elastic peak. 

The measurements with MARI at temperatures be- 
tween 6.2 and 12 K were performed with the c-axis of 
the sample either perpendicular or parallel to the incom- 
ing beam. The dynamic structure factor S{Q,E) was 
measured for wave-vector transfers along the c-axis be- 
tween I 2± 0.175 and ~ 1.5, where the wave- vector along 
the chain is — . Because there are two Ni atoms in each 

c 

unit cell along c, / = 1 corresponds to the antiferromag- 
netic (AF) wave-vector, tt, generally used in theoretical 
work on AF chains. The data with I > 1.5 was contam- 
inated by phonon scattering and so no analysis of the 
magnetic contribution was then attempted. 

FIG. 1. The neutron scattering intensity observed using 
MARI as a function of energy transfer and wave- vector trans- 
fer along the [001] direction. The data shown is an average 
of measurements at 6.2 and 12 K and is given on a loga- 
rithmic scale. The intense band of scattering arises from the 
well-defined Haldane modes and the scattering continuum is 
at higher energies up to ^ 12 meV for 0.6 < I < 1.4. The solid 
lines are the boundary of the continuum that would be ob- 
served for an S—1/2 linear chain if the maximum of the lower 
boundary was at the same energy as the Haldane excitation 
for / = 0.5. 

Constant-Z scans were constructed from the measured 
S{Q,E) by mapping on (/,i?)-space (Fig. p. This data 
was then averaged in stripes with a width — 0.05—0.10 
and corrected for neutron absorption. Fig. ^ shows di- 
rectly the continuum scattering. It is stronger at AF 
momenta around I = 1 than in the region near / = 0.2, 
which is known to carry vanishing magnetic weight pT| , 
and the zone boundary region / = 0.5 A detailed de- 
scription of the analysis of the experiment will be given 
elsewhere Here, we focus only on the scattering 

close to the AF point as measured in the central detec- 
tor bank. As shown in Fig. |^a) the response at Qc = tt 
consists of the sharp mode observed previously but also 
has a component that extends to high frequencies. This 
continuous spectrum has not been observed before. It is 
well above background and moreover is seen to be much 
larger than the scattering near Qc = where the mag- 
netic scattering is expected [|l^ to be vanishingly small. 
We will show below that the continuum cannot arise from 
resolution broadening of the well-defined mode. 

The magnetic excitations in the MARI experiment 
were best fitted by an antisymmetrized Lorentzian peak 
weighted with the detailed balance factor and convoluted 
with the line-shape of the quasi-elastic incoherent scat- 
tering, which is a good estimate of the resolution line- 
shape particularly at low energy transfers. An intrinsic 
Lorentzian peak is supported by a recent theoretical pre- 
diction 1 12 1 and a recent experiment ||l^]. It gives an 
excellent fit to the low energy resonant Haldane peak. 



Nevertheless we see in Fig. that the scattering in- 
tensity at energies above the gap energy is significantly 
higher than expected from the fits to the sharp peak. The 
continuum scattering is clearly observed for 0.6 < I < 1.4 
(Fig. |l|). At 6.2 K, the integrated intensity of the con- 
tinuum is 10(4)% of the total magnetic scattering at the 
AF point. The continuum was also observed in an ex- 
periment in a different configuration at 12 K, where the 
c-axis was parallel to the incoming beam. 

Calculations were made using an ISIS spectrometer 
simulation programme which takes into account the spec- 
trometer parameters such as the detailed pulse line- 
shape, chopper characteristics etc. and predicts the scat- 
tering line-shape for a particular sample orientation and 
model of the scattering cross-section. We find that the 
resolution function reproduced the measured line-shape 
of the quasi-elastic incoherent scattering and of the well- 
defined excitations, but that it cannot give rise to a high- 
energy tail in the region of the continuum. 

FIG. 2. a) Neutron scattering intensity at 6.2 K as a func- 
tion of energy transfer at Z = and I — 1 as measured us- 
ing the time-of-flight spectrometer MARI at ISIS and after 
correction for neutron absorption. The wave-vector trans- 
fer describing the in-plane wave-vector [hhO] varies between 
h — 0.09 and 0.81. The solid line is the fit described in the 
text, b) Neutron scattering counts at 8.5 K as a function of 
energy transfer as measured using the triple-axis spectrome- 
ter DUALSPEC and corrected for neutron absorption. The 
energy scan at the zone boundary (0.31 0.31 1.5) shows that 
the magnetic excitation peak is at ~ 6 meV and that the 
background at low and high energy transfer is considerably 
lower than the neutron scattering at (0.31 0.31 1) between 4 
and 10 meV. The solid lines are fits described in the text in- 
cluding the resolution function. The shaded areas are a guide 
to the eye. 

The continuum scattering was investigated further us- 
ing the triple-axis spectrometer DUALSPEC at Chalk 
River Laboratories and constant-Q scans were performed 
at 8.5 K for various wave- vectors in reciprocal space. The 
energy scans close to the AF point / = 1 revealed consid- 
erable continuum scattering above the well-defined Hal- 
dane excitation. The continuum extends up to 12 meV 
as shown in Fig. ^) for the energy scan at (0.31 0.31 1). 
This wave-vector was chosen to be close to the 3D or- 
dering wave- vector such that the Haldane excitation was 
at the lowest energy to give the largest possible energy 
window for the observation of the continuum. The scat- 
tering intensity at the AF zone boundary, I — 1.5, be- 
low and above the sharp excitation is considerably lower 
than the intensity id I — 1. The magnetic excitation was 
then fitted by a dispersion relation (an antisymmetrized 
Lorentzian weighted by the Bose factor) convoluted with 
the resolution ellipsoid given by Cooper-Nathans's ex- 
pression fl^ j. The dispersion relation has been given by 
and at 6 K the zone boundary energy is 6 meV, the 
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gap at (0.81 0.81 1) is 1.23 meV l|l3| and the bandwidth 
of the gap along [110] direction is given by the gap at 
(0.33 0.33 1), which is 0.35 meV f^]. The Lorentzian 
half width of the excitation is 0.35 meV ||l^]. The fit 
gives an excellent account of the right-hand-side of the 
peak, as shown in Fig. |^), but it cannot explain the 
slowly decreasing scattering continuum at energies above 
4 meV. 

The integrated intensity for the DUALSPEC data at 
the AF point I = 1 was inferred by scaling the peak 
intensity for I — 1.5 according to the /-dependence of 
the intensity measured by the MARI experiment. The 
integrated intensity of the continuum at 8.5 K is found 
to be about 12(3)% of the overall scattering at the AF 
wave- vector. This confirms the results obtained using the 
MARI spectrometer. 

The temperature dependence of the continuum scat- 
tering was further investigated using the RITA triple- 
axis spectrometer at Ris0 and constant-Q scans were 
performed at the ID point (0.81 0.81 1). In Fig. | we 
show that the peak broadens and increases in energy with 
increasing temperature as reported elsewhere p^ . The 
well-defined peak was well fitted by an antisymmetrized 
Lorentzian weighted by the Bose factor and convoluted 
with the resolution ellipsoid. The flat background was 
determined from the scattered intensity below 25 K at 
the highest energy transfers measured {E > 12 meV), 
where the scattering was temperature independent and 
assumed to be non-magnetic. At low temperatures and 
energies higher than the well-defined peak energy, there 
was considerable continuum scattering (Fig. ^). After 
correction for neutron absorption, the total integrated 
intensity at 6 K and 12 K was about 14(3)% and 13(6)% 
of the integrated scattering at / = 1, respectively. At 
higher temperatures, the Lorentzian peak increases in en- 
ergy and broadens so that the continuum could not be 
clearly identified. Nevertheless there was more scattering 
around 8 — 12 meV energy transfer at 6 and 9 K than at 
30 K (see Fig. |). 

FIG. 3. Neutron scattering intensity at / = 1 wave-vector 
transfer as a function of energy transfer measured using the 
triple-axis spectrometer RITA in Ris0 for three different tem- 
peratures. The data is corrected for neutron absorption. The 
solid lines are the fits explained in the text including the spec- 
trometer resolution and the shaded areas are a guide to the 
eye. 

There are a number of theoretical predictions for the 
strength of a continuum at Z = 1. A numerical diagonal- 
ization of a spin-1 chain with iV = 20 sites and nearest- 
neighbor exchange interaction showed that at zero tem- 
perature the continuum above the Haldane excitation 
carries 3% of the total weight ||l^. In the non-linear 
sigma model (NLctM), a continuum arises from 3-particle 
scattering, and its integrated intensity is about 2% of the 



spectral weight ||T^,|l8| . If the coupling of the magnetic 
chain is taken into account via a Random Phase Approx- 
imation (RPA) , the 3-particle scattering becomes depen- 
dent on the wave- vector transfer perpendicular to the ID 
axis, but its integrated intensity remains between 0.5 and 
2.25% of the total spectral weight |l|]. 

The portion of the intensity in the continuum scatter- 
ing for CsNiCls is about 12(2)% for I = 1. Our measure- 
ments for CsNiCla show that the intensity between I = 1 
and I = 1.5 reduces by a factor of 5.4, which is substan- 
tially lower than the theoretical prediction of 8.9. The 
absolute intensity of continuum scattering may therefore 
be only 7(2)% of the theoretical predicted intensity, but 
this is nevertheless much greater than the theoretical pre- 
dictions. 

A considerably stronger 3-particle continuum of about 
17% is predicted by the Majorana Fermion Theory 
(MFT), which describes ID magnetic systems in terms 
of a perturbation theory about the Heisenberg point of 
a model with biquadratic interactions as strong as the 
bilinear exchange jl^. It is unclear whether this model 
is applicable for CsNiCls. Furthermore, no indication of 
the pronounced peak at four times the gap energy, which 
is predicted for this model, was found in our measure- 
ments. Finally the MFT results are not in agreement 
with the numerical calculations of Takahashi |l^ . 

In order to confirm that the continuum scattering 
arises from 3-particle scattering as predicted by the 
NLcrM we have searched for the onset of the continuum at 
3A and for the pronounced maximum of the continuum 
at 6A as predicted by the NLcrM. With A = 0.94 meV, 
our resolution is more than adequate to observe an in- 
creasing continuum intensity between 3 A — 2.8 meV 
and 6A = 5.6 meV, but instead the spectrum shows a 
steadily decreasing intensity (Fig. |2[ ^). This may arise 
because the 3D interactions in CsNiCls cause a change in 
the excitation spectrum which is not properly taken into 
account via RPA and the pronounced maximum at 6 A 
gets much broader as expected from kinematics. Possibly 
the effect can also account for the observed continuum 
intensity being higher than predicted for a ID chain. 

The observed continuum. Fig. |l|, ||, |[ has intensity 
extending to approximately 12 meV, which is twice the 
maximum 1-particle energy for CsNiCla . Because similar 
measurements near / = 0, 2 show no evidence of scatter- 
ing, we conclude that the continuum is strongest in a 
broad region where the AF fluctuations are largest. It 
is then possibly worth commenting that qualitatively the 
continuum scattering has features that are similar to the 
continuum found for 5=1/2 systems (Fig. |^) 

It is clear that further theoretical work is needed to 
understand these results. It is unclear flrstly whether a 
quantitative theory can be obtained from the 3-particle 
picture of the NLcrM or whether it is more appropriate to 
use another development such as the MFT approximation 
pl[ to the Heisenberg model with biquadratic exchange 
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and secondly whether the form of the scattering depends 
on the inter-chain couphng or is characteristic of coupled 
chains. 

In summary, we have measured the continuum scat- 
tering of the spin-1 chain compound CsNiCla in a broad 
region that favours antiferromagnetic fluctuations and 
shown that it is stronger than expected for an ideal spin-1 
chain both from numerical diagonalization and field the- 
oretical calculations involving the NLcrM. These results 
are presented as observations from experiment, and while 
we cannot explain them we hope that they will stimulate 
further theoretical work. 
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